Continuous theory of active matter systems with metric-free interactions 



o 



O 



Anton Peshkov, 1,2 Sandrinc Ngo, 1 Eric Bcrtin, 3 Hugucs Chatc, 1 and Francesco Ginelli 4,5 

'Service de Physique de I'Etat Condense, CEA-Saclay, 91191 Gif-sur-Yvette, France 

2 LPTMC, CNRS-UMR 7600, Universite Pierre et Marie Curie, 75252 Paris, France 

5 Universite de Lyon, Laboratoire de Physique, ENS Lyon, CNRS, 46 Alice d'ltalie, 69007 Lyon, France 

^Istituto dei Sistemi Complessi, CNR, via dei Taurini 19, 1-00185 Roma, Italy 

5 Institute for Complex Systems and Mathematical Biology, King's College, 

University of Aberdeen, Aberdeen AB24 SUE, United Kingdom 

(Dated: March 1, 2013) 

We derive a hydrodynamic description of metric-free active matter: starting from self-propelled 
particles aligning with neighbors defined by "topological" rules, not metric zones, — a situation 
advocated recently to be relevant for bird flocks, fish schools, and crowds — we use a kinetic approach 
to obtain well-controlled nonlinear field equations. We show that the density-independent collision 
rate per particle characteristic of topological interactions suppresses the linear instability of the 
homogeneous ordered phase and the nonlinear density segregation generically present near threshold 
in metric models, in agreement with microscopic simulations. 
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Collective motion is a central theme in the emerging 
field of active matter studies [1(. For physicists, the inter- 
est largely lies in the non-trivial cases where the emer- 
gence of collective motion can be seen as an instance 
of spontaneous symmetry-breaking out of equilibrium: 
without leaders, guiding external fields, or confinement 
by boundaries, large groups inside which an "individual" 
can only perceive local neighbors are able to move co- 
herently After this was realized in the seminal papers of 
Vicsek et al. [2| and Toner and Tu |3( , much progress has 
been recorded in the physics community [l|, |j] , alongside 
continuing modeling work in ethology and biology [a, lfj| . 

Most models consist of self-propelled particles interact- 
ing with neighbors defined to be those particles within 
some finite distance [7J. Among those "metric mod- 
els", that introduced by Vicsek et al. [2[ is arguably 
the simplest: in the presence of noise, point particles 
move at constant speed, aligning ferromagnetically with 
others currently within unit distance. The study of 
the Vicsek model has revealed rather unexpected be- 
havior. Of particular importance in the following is the 
emergence of phase segregation, under the form of high- 
density high-order traveling bands [8|, in a large part 
of the orientationally-ordercd phase bordering the on- 
set of collective motion, relegating further the spatially- 
homogeneous fluctuating phase treated by Toner and Tu. 
Similar observations of density segregation were made for 
important variants of the Vicsek model, such as polar 
particles with nematic alignment JfJ (self-propelled rods) 



or the active nematics model of [lfj, lll| . The genericity 
of these observations has been confirmed, in the Vicsek 
case, by the derivation and analysis of continuous field 
equations 12|,|l3| (see also [lj, [l5 1 ) . It was shown in par- 
ticular that the homogeneous ordered solution is linearly 
unstable near onset, and that solitary wave structures 
akin to the traveling bands, arise at the nonlinear level. 



Even though metric interaction zones are certainly of 
value in cases such as shaken granular media [lot Il7| and 
motility assays [lg, [l9( where alignment arises mostly 
from inelastic collisions, it has been argued recently |20l — 
[22| that they are not realistic in the context of higher 
organisms such as birds, fish, or pedestrians, whose nav- 
igation decisions are likely to rely on interactions with 
neighbors defined using metric-free, "topological" crite- 
ria. Statistical analysis of flocks of hundreds to a few 
thousand individuals revealed that a typical starling in- 
teracts mostly with its 7 or 8 closest neighbors, regardless 
of the flock density [20( . The realistic, data-based, model 
of pedestrian motion developed by Moussaid et al. relies 
on the "angular perception landscape" formed by neigh- 
bors screening out others |2lj . 

At a more theoretical level, the study of the Vicsek 
model with Voronoi neighbors [23| (those whose associ- 
ated Voronoi cells form the first shell around the cen- 
tral cell) has shown that metric-free interactions are rel- 
evant at the collective scale: in particular, the traveling 
bands mentioned above disappear, leaving only a Toner- 
Tu-like phase. Below, we show that the introduction 
of Voronoi neighbors suppresses the density-segregated 
phase in other variants of the Vicsek model. In spite of 
the recognized importance of metric-free interactions, no 
continuous field equations describing the above models 
are available which would help put the above findings on 
firmer theoretical ground. 

In this Letter, starting from Vicsek-style microscopic 
models with Voronoi neighbors, we derive nonlinear field 
equations for active matter with metric-free interactions 
using a kinetic approach well-controlled near the on- 
set of orientational order. We show that the density- 
independent collision rate per particle characteristic of 
these systems suppresses the linear instability of the ho- 
mogeneous ordered phase and the nonlinear density seg- 



rogation in agreement with microscopic simulations. We 
finally discuss the consequences of our findings for the 
relevance of metric-free interactions. 

Let us first stress that with metric-free interactions, 
say with Voronoi neighbors, the tenet of the Boltzmann 
equation approach [24[ , the assumption that the system is 
dilute enough so that it is dominated by binary collisions, 
is never justified since, after all, a particle is almost con- 
stantly interacting with the same number of neighbors. 
Here, instead, we introduce an interaction rate per unit 
time. In the low rate limit, binary interactions dominate, 
and one can proceed "as usual" . 

Our starting point is not a Vicsek-style model, but a 
microscopic rule in the same universality class: N point 
particles move at constant speed tj on a L x L torus; 
their heading is submitted to two different dynamical 
mechanisms, "self-diffusion" and aligning binary "colli- 
sions" . In self-diffusion, 8 is changed into 8' = 8 + r\ 
with a probability A per unit time, where 77 is a ran- 
dom variable drawn from a symmetric distribution P„ (77) 
of variance a 2 . Aligning "collisions" occur at rate a 
per unit time with a randomly chosen Voronoi neighbor. 
Then 9 is changed to 8' = *(0, n ) + 77 where n is the 
heading of the chosen neighbor and the noise 77 is also 
drawn, for simplicity, from P a (n). Isotropy is assumed, 
namely *(0i + 0, 8 2 + <j>) = *(0i, 62) + <j> [2ir]. For the 
case of ferromagnetic alignment treated in detail below, 
^(0i,0 2 ) =arg(e iei W 2 ). 

The evolution of the one-particle phase-space distribu- 
tion f(r,8,t) (defined over some suitable coarse-grained 
scales) is governed by the Boltzmann equation 

d t f(r, 0, t) + v e(0) ■ V/(r, 0, t) = I diS [/] + 7 co ii[/] (1) 

where e(0) is the unit vector along 8. The self-diffusion 
integral is 

/■7T /"OO 

Jdiff[/] = -A/(0) + A/ old' di 1 P a (i 1 )S 27r (8'-8 + r 1 )f(9') 
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(2) 

where S 27T is a generalized Dirac delta imposing that the 
argument is equal to zero modulo 2tt. In the small- 
a limit, orientations are decorrelated between collisions 
("molecular chaos hypothesis"), and one can write: 

/7T /*7T POO 

dej dd 2 dr, 

x P CT (T7)/(0i)/(02)M*(0i>02)-0+r?) (3) 

The main difference with the metric case treated in [12l | is 
the "collision kernel" , which is independent from relative 
angles and inversely proportional to the local density 



p(v,t) = / f(r,9,t)M. 



(4) 



Note that, in agreement with basic properties of mod- 
els with metric- free interactions, Eq. ([!}, together with 



definitions ([2]|4]), is left unchanged by an arbitrary nor- 
malization of / (and thus of p) and thus does not de- 
pend on the global density po = N/L 2 . Furthermore, a 
rescaling of time and space allows to set A = vq = 1, a 
normalization we adopt in the following. 

Equations for hydrodynamic fields are obtained by ex- 
panding /(r, 8,t) in Fourier series, yielding the Fourier 
modes f k (x,t) = J* 7T d9f(x,9,t)e ike where f k and f_ k 

arc complex conjugates, /o = p. and the real and imag- 
inary parts of /1 are the coordinates of the momentum 
vector w = pP with P the polar order parameter field. 
Using these Fourier modes, the Boltzmann equation ([T]) 
yields an infinite hierarchy: 
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where the complex operators V = d x 
d x — idy have been used, the binary collision rate a is now 
expressed in the rescaled units, P k = J^ c co d7jP a (n)e' lkri 
is the Fourier transform of P a , and J kq is an integral 
depending on the alignment rule *5f. Below, we specialize 
to the case of ferromagnetic alignment, for which 

For k — the r.h.s. of Eq. ([5]) vanishes and one recovers 
the continuity equation: 



d t p + V • w = 0. 



(7) 



To truncate and close this hierarchy, we assume the fol- 
lowing scaling structure, valid near onset of polar order, 
assuming, in a Ginzburg-Landau-likc approach, small 
and slow variations of fields: 
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Note that the scaling of space and time is in line with 
the propagative structure of our system [25( . The lowest 
order yielding non-trivial, well-behaved equations is e 3 : 
keeping only terms up to this order, equations for f k - >2 
identically vanish, while f 2 , being slaved to /1, allows 
to close the /j equation which reads, in terms of the 
momentum field w: 



<9 t w 



7(w • V)w = — — Vp ■ 
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Apart from some higher order terms wc have discarded 
here, this equation has the same form as the one derived 
in [121 ] for metric interactions, but with different trans- 
port coefficients: 
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Note first that, contrary to the metric case, the coefficient 
p of the linear term does not depend on the local density 
p; coefficients of the non-linear terms depend on density 
to compensate the density dependence of w. Note fur- 
ther that v, k, and £ are positive since < Pfc < 1, 
so that in particular the nonlinear cubic term is always 
stabilizing. For an easier discussion, we consider now 
the Gaussian distribution P a (il) = — ^= ex p[ — ^s] f° r 

which Pk — exp[— k 2 a 2 /2}. Then p is negative for large 
a (where the trivial w = solution is stable with respect 
to linear perturbations), and changes sign for [26| 
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For <t < er c , the nontrivial homogeneous solution p = po, 
w = wi = e vm/£ (where e is an arbitrary unit vector) 
exists and is stable to homogeneous perturbations. 

We now focus on the linear stability of wi with respect 
to arbitrary wavevector q. Because we want to discuss 
later differences between the metric and metric-free cases, 
we keep a formal /^-dependence of the linear transport 
coefficient. Linearization around wi yields 



dtSp = — V • (5w 



(12) 



d t Sw = — 7(w 1 -V)<5w — V(5p+^V 2 i5w-|-«;V(w 1 -i5w) 
— kwi(V-<5w) — 2£wi(wi-<5w) + (// — £'w 2 ) wiSp 

where primes indicate derivation with respect to p. Using 
the ansatz (<5/o(x, i), <Jw(x, t)) = exp(si + iq-x)(<5p q ,<5w q ) 
allows to recast Eq. (fT2j) as an eigenvalue problem for s. 

Wc have solved numerically this cubic problem for the 
metric-free case using coefficients (|10[) and Gaussian noise 
in the full (a, a) parameter plane. The resulting stabil- 
ity diagram, presented in Fig. [TJ shows that, contrary 
to the metric case, the homogeneous ordered phase is 
stable near onset. Like in the metric case 27|, there ex- 
ists an instability region to oblique wavevectors of large 
modulus rather far from the transition line. Given that 
microscopic simulations show no sign of similar instabili- 
ties, we believe that the existence of this region, situated 
away from the validity domain of our approximations, is 
an artifact of our truncation ansatz. 

At the nonlinear level, we performed numerical sim- 
ulations [28| of Eq. ([9]) (again with coefficients (fT0|) and 
Gaussian noise) starting from initial conditions with large 
variations of both p and w. With parameters a and a 2 
in the ordered stable region of Fig.Q] we always observed 
relaxation towards the linearly-stable homogeneous solu- 
tion wi, albeit after typically long transients. Starting in 
the unstable region, the solution blows up in finite-time, 
signaling that indeed our equation is ill-behaved when 
considered too far away from onset. 

The stabilization of the near-threshold region by 
metric-free interactions can be directly traced back to 
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FIG. 1: (color online) Continuous theory for self-propelled 
particles aligning ferromagnetically (Eqs. J7} and ((9j with co- 
efficients (|10|) . Gaussian noise), (a) Phase diagram in the 
(a, a 2 ) plane. The solid line marks the order-disorder tran- 
sition. The homogeneous ordered solution wi exists below 
this line and is linearly stable above the coloured linear in- 
stability region. The colour scale (in radians) indicates the 
most unstable wavevector direction cj>. (b) Modulus q of the 
most unstable wavevector (green line) and the real part of its 
corresponding eigenvalue s+ (dashed red line) as a function 
of a 2 at a = 2. 



the absence of p dependence of p, in agreement with re- 
|29| where the long-wavelength instability 
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of wi was linked to p' > 0. In the long wavelength limit 
q = |q| <C 1, the eigenvalue problem can be solved an- 
alytically with relative ease. The growth rate s is the 
solution of the cubic equation 

s 3 + f3 2 s 2 + /3 lS + A, = (13) 

where the coefficients, to lowest orders in q, are given by 
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and <j) is the angle between q and Wi (which has been 
chosen parallel to the abscissa). Near threshold, where 
our truncation is legitimate and p ~ e 2 , two eigenvalues 
are always stable and linear stability is controlled by the 
dominant solution real part 
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with h(0) = 2 and h(<j>) = 1 otherwise |30(. This expres- 
sion immediately shows that s+ < in the metric-free 
case where pi = 0, confirming the stability of the homo- 
geneous ordered solution wi. This stabilizing effect can 



be ultimately traced back to the (negative) pressure term 
—Vp appearing in Eq. ([9]). Conversely, in the metric case 
for which /i ~ + near threshold, s + is always positive, 
yielding the generic long wavelength instability leading 
to density segregation in metric models. 

To sum up at this point: Through a kinetic approach, 
we have derived a continuous, "hydrodynamic" theory for 
self-propelled particles aligning with topological neigh- 
bors ( "metric- free interactions" ) . Calculations were pre- 
sented in full for the case of ferromagnetic alignment, for 
which we showed that the generic long-wavelength in- 
stability of the homogeneous ordered solution present in 
metric models near onset is suppressed. At the nonlin- 
ear level, we observed numerically that the density seg- 
regated solutions of metric models vanish, so that the 
homogenous solution seems to be a global attractor. 

These results rest on the independence of the linear 
coefficient /i on p, a direct consequence of the fact that 
the interaction rate per particle in topological models is 
fixed only by geometrical constraints and does not grow 
with local density. This property actually holds for any 
metric-free system: all the linear coefficients fi^ appear- 
ing in the equations for ff. read 



/i fc = P fe - 1 - a + aP k (J kk + J, 



ko) 



(16) 



and arc thus independent of p. 

Our other general conclusions also extend to the other 
basic classes of simple, Vicsek-style collective motion: we 
have in particular worked out the case of polar parti- 
cles with nematic alignment ("self-propelled rods") for 
which #(0i, 2 ) = arg(e jei +sign[cos(6q-(9 2 )]e^ 2 ). While 
full details will be given in [3l|, we only sketch here 
the salient points. Let us first recall that with nematic 
alignment, the metric model studied at the microscopic 
level in [9(, shows global nematic order. In a large re- 
gion of parameter space bordering onset, order is segre- 
gated to a high-density stationary band oriented along 
it. We have studied numerically the metric-free version 
of that Vicsek-style model with Voronoi neighbors. Like 
for its ferromagnetic counterpart, no segregation in bands 
is observed anymore, and the transition to nematic order 
is then continuous (Fig. [2]). These properties are well- 
captured, both in the metric and metric-free case, by a 
controlled hydrodynamic approach of the type presented 
here [3l[ . The nematic symmetry of the problem requires 
to consider three hydrodynamic fields 15| , corresponding 
to the modes k — 0, 1, 2 in Eq. ([5]), with fa coding for the 
nematic tensor field pQ. We have performed the analy- 
sis of the 5x5 linear problem expressing the stability of 
the homogeneous nematically-ordered solution (w = 0, 
/?Q = Cst.) appearing at onset in both the metric and 
non- metric cases |31[. Whereas the metric case shows 
a long wavelength, transversal instability of the homoge- 
neous ordered solution near onset, this solution is linearly 
stable in the metric-free case. Again, this difference can 
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FIG. 2: (color online) Vicsek-style model with nematic align- 
ment and topological neighbors, where N — L 2 particles 
move at speed vo — i on a L x L torus. Headings and 
positions are updated at discrete timesteps according to 

0] +1 = arg [E fe ^ Sign[cos(0| - e))]e< + n)a$\ and v] +1 = 
Tj +vo e(#' +1 ) where e(6) is the unit vector along 0, the sum 
is over the n* Voronoi neighbors of particle j (including j it- 
self), and £j is a random unit vector in the complex plane. 

Global nematic order parameter S = (j-^ ^2 k e t29k \) t (a) and 
its Binder cumulant G (b) vs a for L — 32, 64, 128 (the arrows 
indicate increasing sizes), revealing a continuous transition. 



be traced back to the p-dependence of the linear coef- 
ficients [i k . At the nonlinear level, simulations indicate 
that the homogeneous ordered solution is a global attrac- 
tor in the metric-free case. 

Our analysis can also be extended to "non-ballistic" 
active matter such as the driven granular rods model 
( "active nematics" ) studied in lfj, [ll( which, for metric 
interactions, also shows near-threshold phase segregation 
[lot |32( . Simulations of the metric-free microscopic ver- 
sion (with Voronoi neighbors) show no such segregation. 
In a kinetic approach, because active nematic particles 
move by non-equilibrium diffusive currents rather than 
by ballistic motion, the Boltzmann equation has to be re- 
placed by a more general master equation. But it is nev- 
ertheless possible to derive a continuous theory which, in 
the metric-free case, yields a homogeneous ordered phase 
stable near onset for essentially the same reasons as in 
the cases presented above [331 ] . 

In conclusion, simple, Vicsek-style, models of active 
matter where self-propelled particles interact with neigh- 
bors defined via non-metric rules (e.g. Voronoi neigh- 
bors) are amenable, like their "metric" counterparts, to 
the construction of continuous hydrodynamic theories 
well-controlled near onset. The relatively simple frame- 
work of Vicsek-style models offers a two-dimensional pa- 
rameter plane which can be studied entirely. More com- 
plicated microscopic starting points, for instance where 
positional diffusion would also be considered, inevitably 
raise the dimensionality of parameter space. We have 
shown here that the non-metric theories differ essentially 
from the metric ones in the independence of their linear 
coefficients fi k on the local density, a property directly 
linked to the fact that the collision rate per particle is 



constant in metric-free systems. We have shown further 
that the homogeneous ordered phase is linearly stable 
near onset for metric-free systems, in contrast with the 
long-wavelength instability present in metric cases. 

We finally discuss the relevance (say in the 
renormalization-group sense) of metric-free interactions 
in deciding active matter universality classes. Our 
work has shown that the deterministic continuous the- 
ories of metric-free active matter systems are formally 
equivalent to those of their metric counterparts, ex- 
cept for the density-dependence of the linear coefficients. 
This could be taken as an indication, in the case of 
ferromagnetically- aligning particles, that the homoge- 
neous, ordered, fluctuating phase observed in the Vic- 
sek model with Voronoi neighbors docs not differ from 
the Toner- Tu phase of its metric counterpart, in con- 
tradiction with the numerical discrepancies between the 
two cases reported in 23[ about the scaling exponent 
of the anomalously-strong density fluctuations. This 
calls for more extensive microscopic simulations assess- 
ing finite-size effects, but also for incorporating effective 
noise terms, properly-derived in both cases, and to study 
the resulting field theories, a task left for future work. 
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Marchctti, J. Prost, and J. Toner for enlightening discus- 
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Group 2011/2012: Statistical Physics of Collective Mo- 
tion. 
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